c_big=c_c_big;
bp=bp_cc;
TempReshape=repmat(c_def_2d(:),1,NBG)'; 
c_def=reshape(TempReshape(:),NBG,NB,NSS); % THIS SHAPES TO GET THE FINAL THING
c_big_bad=c_def_2d;
bp_bad=bp_def_2d;

 parfor j=1:NSS %Today's Exogenous Shock
           for t=1:NB %Today's private debt index 
                for i=1:NBG %Today's public Debt shock
                        price_def(i,t,j)=((1-omega)/omega*(c_def(i,t,j)/yN(j))^(1+ita));
                        for inext=1:NBG %This is the Grid Search                            
                            price_cc(i,t,j,inext)=((1-omega)/omega*(c_c_big(i,t,j,inext)/yN(j))^(1+ita));
                        end 
                 end
           end
 end